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Abstract 

In this note we present numerical simulations of binary mixtures and we find indica- 
tions for a thermodynamic transition to a glassy phase. We find that below the transition 
point the off equilibrium correlation functions and response functions seems to be asymp- 
totically compatible with the relations that were originally derived Cugliandolo Kurchan 
for generalized spin glasses. 

1 Introduction 

The behaviour of an Hamiltonian system (with dissipative dynamics) approaching equilibrium 
is well understood in a mean field approach for infinite range disordered systems |], §, |3| . In this 
case we must distinguish an high and low temperature region. In the low temperature phase 
the correlation and response functions satisfy some simple relations derived by Cugliandolo 
e Kurchan [|IJ. We present in this note the first investigation of the relations among these 
quantities for binary glasses. We find indications for the existence of a phase transition. The 
numerical data are compared to the theory. The results of this comparison point toward the 
applicability of the results of the CK dynamical theory to binary glasses. 

Generally speaking in a non equilibrium system it is natural to investigate the properties 
of the correlation functions and of the response function. Let us concentrate our attention 
on a quantity A(t), which depends on the dynamical variables x(t). Later on we will make a 
precise choice of the function A. 

Let us suppose that the system starts at time t — from a given initial condition and 
subsequently it follow the laws of the evolution at a given temperature T. If the initial 
configuration is not at equilibrium at the temperature T, the system will display an off- 
equilibrium behaviour. In many case the initial configuration is at equilibrium at a temperature 
X" > T; different results will be obtained as function of T' . In this note we will consider only 
the case T' » T (in particular we will study the case T' = oo). 

We can define a correlation function 

C{t,t w )=<A{t w )A{t + t w )> (1) 
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and the response function 

5A(t + 1 % 



G(t, t w ) = 



e=0 



Se(t t 

where we are considering the evolution in presence of a time dependent Hamiltonian in which 
we have added the term 

' dte(t)A(t). (3) 



The off-equilibrium fluctuation dissipation theorem of the CK theory states some properties 
of the correlation functions in the limit t w going to infinity. The usual equilibrium fluctuation 
dissipation (FDT) relation tell us that 



where 



G(t)= lim G(t,t w ), C(t)= lim C(t,t w ). (5) 

tin — KX> tin — >oo 



In our notation the correlation and response functions which depend on two times are the 
off-equilibrium ones; those which depend on only one time are the equilibrium ones, i.e. they 
are the limit of the off-equilibrium correlations when both times goes to at infinity at fixed 
distance. 

It is convenient to define the integrated response: 



R{t,t w ) = f drG(r,t, 
Jo 



R(t) = lim R(t,t w ), (6) 

t-u, — >0O 

which is the response of the system to a field acting for a time t. 
We can also define the quantities 

C(oo) = lim C(t), R(oo) = lim R(t). (7) 

t — >oo t — ►oo 

The static FDT relation is 

R(oo) =P(C(0) -C(oo)). (8) 
The last equation can be naively written also as 

R(oo) = (3 (< A 2 > - < A > 2 ) , (9) 

Here the brakets denote the usual equilibrium expectation value. If there is only one equilib- 
rium state (or two that have opposite values < A > and differ by a symmetry), the previous 
formula, eq. (|9]) , is correct, otherwise a more lengthy discussione is need. 

A very interesting situation happens when the quantity < A > is identically zero because of 
symmetry arguments in the high temperature. It is quite possible that there is a spontaneous 
symmetry breaking: two or more states of the systems may be present and the expectation 
of A in the appropriate state becomes different from zero. A typical example of this situation 
(see the appendix for more details) is given by spin glasses 0, |^, |6|, where the magnetic 
susceptibility can be written at zero magnetic field as x — P/NJ2i=i jv(< °f > — < °i >2 ) = 
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In this case the following relation is valid in the high temperature phase 

R(oa)=pC(0), (10) 

where we recall that in our notations C(0) = lim^oo C(t,t). 

The breaking at low temperature of the relation eq. (|T0"| ) is a signal of a phase transition. 
We will denote by T c the temperature at which the previous relation breaks. We can also 
introduce an order parameter defined by 

— 1 P^D n 1 x 

R[oo) 

where 

C D = lim C(t,t). (12) 

t — >oo 

We can get further information on the nature of the transition if we stay in the framework of 
the CK theory for the approach to equilibrium . In the study of off-equilibrium spin glasses 
systems Cugliandolo and Kurchan proposed that the response function and the correlation 
function satisfy the following relations: 

G(t, t w ) « -0X (C(t, t w )) ggfeM , (13) 

ot 

The previous relation can be also written in the following form 

R(t,t w ) « (3 / X{C)dC. (14) 

Jc(t,t w ) 

The function X(C) is system dependent and its form tell us many interesting information. 
If C(oo) 7^ 0, we must distinguish two regions: 

• A short time region where X = 1 (the so called FDT region) and C > C(oo). 

• A large time region (usually t = 0{t w ) where C < C(oo) and X < 1 (the aging region) 

EE- 

In the simplest non trivial case, i.e. one step replica symmetry breaking, the function 
X{C) is piecewise constant, i.e. 

X{C) = m for C < C(oo), 
X{C) = 1 for C > C(oo). (15) 

In all known cases in which one step replica symmetry holds, the quantity m vanishes 
linear with the temperature at small temperature. It often happens (but it is not compulsory) 
that m = 1 at T = T c . 

We notice that we must be quite careful in exchange limits in the low temperature phase: 
the correlation function C satisfy the relation 

lim C(t,t w ) = ^ lim ( lim C(t,t w )) = C(oo). (16) 

t—>-oo t— »oo \t w — >oo / 
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In the sane way we have that in the region where t and t w are both large 
R(t, t w ) « R(oo) = /3(C(0) - C(oo)) for t « ^, 

R(t,t w ) ^ R eq = (3 [° dCX(C) for t»t w . (17) 

Therefore it is quite possible that 

lim R(t, t w ) = R eq Moo). (18) 

This phenomenon is present as soon the function X(C), is non zero outside the FDT region 
and it is the typical situation that happens when replica symmetry is broken [[|, [| . 

The previous considerations are quite general. However the function X(C) is system 
dependent and its form tell us many interesting information. Systems in which the replica 
symmetry is not broken are characterized by having m = in the formula eq. (fL5 ) . 



Sometime simple aging is also assumed f7j, i.e. the following the scaling relation holds 
outside the FDT region: 

C(t,t w ) = C s (pj. (19) 

Simple aging may be correct, but it is not a necessary consequence of the previous relations 
and its verification is not the primary aim of this note (a discussion of simple aging in the 
same system can be found in M). 

The aim of this note is to show that binary mixture of spheres do have a transition in the 
thermodynamics sense at a temperature near the glassy transition, which can be characterized 
by a non zero value of the appropriate order parameter q^. More precisely we will show that 
there are equilibrium quantities which have an irregular (i.e. non-analytic behaviour) at the 
transition point T c . Moreover the correlation and response functions seem to satisfy the 
relations of the CK theory, with the function X(C) is compatible to be given by the one step 



formula (15) 



The paper is organized as follows. In section II, we define the model, the relevant quantities 
(i.e. the asymmetry or stress) and we present some general considerations. In section III we 
study the approach to equilibrium of quantities defined at given time, e.g the energy and the 
equal time fluctuations of the stress. We show that the fluctuations of the stress are strongly 
indicative of a phase transition. In section IV we make a comparison of our data with the CK 
theory of aging. Finally in the last section we present our conclusion. At the end of the paper 
there is a short appendix where some results on spin glasses are recalled and a comparison is 
done with the finding of this paper. 



2 The model 

2.1 The Hamiltonian 

The model we consider is the following. We have taken a mixture of soft particles of different 
sizes. Half of the particles are of type A, half of type B and the interaction among the particle 
is given by the Hamiltonian: 
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where the radius (a) depends on the type of particles. This model has been carefully studied in 
the past ]n| frT[ [12], p~3|| . It is known that a choice of the radius such that o B jo A = 1.2 strongly 
inhibits crystallisation and the systems goes into a glassy phase when it is cooled. Using the 
same conventions of the previous investigators we consider particles of average diameter 1, 
more precisely we set 

a\ + 2(a A + a B f + a% = ^ ^ 

Due to the simple scaling behaviour of the potential, the thermodynamic quantities depend 
only on the quantity T 4 / p, T and p being respectively the temperature and the density. For 
defmiteness we have taken p — 1. The model as been widely studied especially for this choice 
of the parameters. It is usual to introduce the quantity r = /5 4 . The glass transition is known 



to happen around V = 1.45 [11 



Our simulation are done using a Monte Carlo algorithm, which is more easy to deal with 
than molecular dynamics, if we change the temperature in an abrupt way. Each particle is 
shifted by a random amount at each step, and the size of the shift is fixed by the condition that 
the average acceptance rate of the proposal change is about .4. Particles are placed in a cubic 
box with periodic boundary conditions. In our simulations we have considered a relatively 
small number of particles N = 18, N = 34 and N = 66. We start by placing the particles 
at random and we quench the system by putting it at final temperature (i.e. infinite cooling 
rate). 



2.2 The stress 

The main quantity on which we will concentrate our attention is the asymmetry in the energy 
(or stress): 

A = J2 frf'J + ^y 2 (2(a, - x k f - { Vi - y k f - ( Zi - z k f) = 2T M - T 2 , 2 - T 3 , 3 . (22) 



i<k 



In other words A is a combination of the diagonal components of the stress energy tensor. If 
the particles are in a cubic symmetric box, we have that 

< A >= 0. (23) 

If the box does not have a cubic symmetry the effect of the boundary disappears in the infinite 
volume limit (at fixed shape of the boundary) and we have that the stress density a vanishes 
in this limit: 

< A > 

lim =o = 0. (24) 

iV^oo N 

What happens when we add a term eA to the Hamiltonian is remarkable. Let us consider 
the new Hamiltonian 

H + !2eA, (25) 

where H was the old Hamiltonian eq. (p0|) . 

It is convenient to consider the following Hamiltonian: 

H ^im±m\ (26) 



i<k 
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where 

r 4 , fc (e) 2 = (x< - x k ) 2 (l + e)- 4 + (fo - y k f + (z, - z k f) (1 + e) 2 . (27) 

It is evident that we can recover the original Hamiltonian by contracting (for positive e) the x 
direction by a factor (1 + e)~ 2 and expanding the y and z directions by a factor 1 + e, keeping 
in this way the volume constant. 

As far as the expectation values of intensive quantities do not depend on the shape of the 
box, we can compute the properties of the theory with e 7^ in terms of those at e = 0. For 
example one finds that the energy and stress density are given by 

e(e)=e(0) o(e) = + 0(e 2 ). (28) 

o 

At the order e the Hamiltonian in eq. ( p6|) coincide with the previous one eq. ( 25|) . 
We thus arrive to the following conclusion, if we consider the response of the system to 
adding an asymmetric term in the Hamiltonian. 

• The equilibrium response function R eq is exactly given by ^e(O) at all temperatures. 

• If there is only one equilibrium state, and this happens in the high temperature phase 
we must have at e = 

<A2> V (29) 



If we define (at e = 0) 



N 5(3 



C (M.) = 12/3 < ^ +t ; ) ^ )> . (30) 



(where the factor 12/5 has been added to simplify the FDT relation), we obtain that in 
the high temperature phase 

36 

limC(M) = C D = — e. (31) 

t->oo 5 

It is convenient to define the quantity 

w ® - W < 32 > 

and investigate its limit for large time (W = lim^oo W(t)). In the same way as in spin glasses 
we can define a quantity by 

q A = l-^ = l-W- 1 (33) 

In the next section we shall see that there is transition from q = in the high temperature 
phase to a non zero value of q at low temperatures. 
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3 On the transition 



3.1 General considerations 

We have done simulations for various values of N ranging from N = 18 to N = 66. For 
N = 18 and iV = 34 we have measured the correlations by using 1000 runs with different 
starting point; for iV = 66 we have only 250 samples. The evolution was done using the 
Monte Carlo method, with an acceptance rate fixed around .4. At the end of each Monte 
Carlo sweep all the particles are shifted of the same vector in order to keep the center of 



mass fixed |p^| . This last step in introduced in order to avoid drifting of the center of mass 
and it would be not necessary in molecular dynamics if we start from a configuration at zero 
total momentum. Most of the quantities that we measure (with only one exception, i.e. the 
quantity Q) to be defined later eq. (|35|) are not affected by this shift. 
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Figure 1: The value of E(t) versus y = 

f(-.T) f or different values of T (1.4, 1.6, 1.8 2.0) at 

N = 18. 

Four observations are in order: 

• We need to do the average over many samples in order to decrease the error on the 
correlation C(t,t). Only one run would practically give no information because for this 
quantity the errors do not go to zero when N goes to infinity. 

• The values of N may look rather small. On the other hand our task is to show the 
existence of a thermodynamic phase transition. Although definite conclusions can be 
only obtained by a careful finite size analysis for large N, strong indications of the 
existence of a transition can be obtained also for small samples. For example in the Ising 
case the study of the susceptibility on a 4 3 lattice is enough to draw strong indications 
for the existence of a transition and the 64 binary degree of freedom are more or less the 
equivalent of our 54 continuous degrees of freedom for N = 18. 
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Figure 2: The value of E(t) versus y = 

£(-•7) f or different values of F (1 .4, 1.6, 1.8 2.0) at 

N = 34. 



The limit t — > oo may be tricky. All the previous theoretical discussion were done for 
an infinite system. In practice we need that that iV is greater that the time to a given 
power, the exponent of the power being system dependent. Therefore we cannot strictly 
take the limit t — > oo for finite N. We have to study the behaviour of the system in 
region of time whose size increase with N. Finite size corrections sometimes become 



much more severe at large times [14 



In this note we have explored a region of not very large times and we have extrapolated 
the data to infinity by assuming simple powers corrections at large times. The real 
situation is likely more complicated. It is quite possible that when we fast cool the 
system we end up in a metastable state of energy higher the equilibrium one and that 
the system is going to relax to the equilibrium values with a much slow process, which 
cannot be seen on this time scale and produces a drift on a much slower scale. In this 
situation our results concern not the real equilibrium value of the energy (and of the 
other quantities), but their value in a metastable state. 



3.2 The energy 

Let us consider the time dependence of the energy. In fig. (HD and (0) we see the energy for 
r = 1.4, 1.6, 1.8, 2.0, which are near or below the value T = 1.44 at which the phase transition 
is estimated from the behaviour of the equilibrium correlation functions. 

The data are plotted as function of y = t^~' 7 '. The data are quite linear as function of y 
at iV = 34, while they show some curvature at large t, small y for N = 18. This seems to be 
a small finite volume effect. Apart from this effect and an overall shift, the two sets of data 
are quite similar so that we have reasons to suppose that the iV = 34 are nearly asymptotic. 
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This is confirmed by a run at N = 66 for T = 1.8: a comparison of the three runs at T = 1.8 
is showed in fig. (|J). 
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Figure 3: The value of E{i) versus y = ' 7 ' for different values of N (18, 34 and 66) at 
r = 1.8. 

The results seem to be puzzling. The energy goes to its asymptotic values with a power 
correction and the exponent of power correction is independent from the temperature when 
the temperature changes of about a factor 4 (and also it is not a small exponent!). On the 
other hand we know that in real systems the convergence of the energy to the asymptotic value 
is extremely slow below the glass transition. It is quite likely that we are blind to this slower 
process that happens on a much slower time scale. Moreover the temperature independence 
of the exponent strongly suggests that the process is not dominated by activated processes 
which become dominant at very larger times. 

The situation is quite reminiscent of the mean field theory case, where the power to ap- 
proach to a metastable state depends weakly on the temperature flT5| . We can thus suppose 
that here also we converge to a metastable state, whose energy may be larger that the equi- 
librium one. It would be interesting to verify this point by carefully computation of the true 
equilibrium energy by more tuned simulations. We can thus tentatively conclude that our 
large time extrapolation do concern some kind of metastable state. If this happens, it is 
quite remarkable that the two time scales (controlling respectively the the approach to the 
metastable state and the slow decay of the metastable state) are so separated that the first 
can be studied independently from the second. 



3.3 The fluctuation of the stress 

Here we perform the same analysis as in the previous subsection for the quantity W(t) defined 
in eq. (|32]), but we are very interested in finding the exact value of extrapolated at infinite 
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time. In the high temperature region we find that W extrapolates to a value very near to 1, 
as expected from the previous considerations. 

New phenomena appear when we decrease the temperature below the one corresponding 
to T = 1.4. In fig. (|) we show the data of W(t) versus y = t (_ - 7) for T = 1.8 for N = 18 and 
N = 66. 
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Figure 4: The value of W(t) versus y = p ' 7 - ) for T = 1.8 for N = 18 and N = 66. 

Here also we the same phenomena as for the energy. On the large sample the data are 
quite linear when plotted versus y. It is remarkable that the same choice of power which (i.e. 
.7) which works for the energy, is also good for W. It seems that we have the same power 
corrections for both quantities. 

The data at small lattice show a curvature in the plot versus y, which is absent in larger 
sample. This is finite volume effect which indicates which is the range of time for which we 
can comfortable assume that the volume is sufficiently large. 

We show in fig. (^) the value of W(t) versus y = t* - 7 ) for different values of T at N = 34. 
A liner fit is rather good and the extrapolated value of W(oo)~ l are shown in fig. (|6]). 

It is clear from both fig. (|5|) and (|5]) that W(oo) becomes different from 1 in the low 
temperature region. If we assume that it diverges at low temperature, the data are compatible 
with the usual situation where W(oo) is roughly proportional to the temperature at small 
temperature. This results should be taken as an indication that, at least in the metastable 
region there is a transition to a region where W(oo) ^ 1. 



4 The comparison with the CK Theory 
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Figure 5: The value of W(t) versus y = t { --^ for different values of V (1.0, 1.2, 1.4, 1.6, 1.6 
and 2) at iV = 34. 
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Figure 6: The value of 1 — qA = ^(oo)™ 1 as function of the temperature using the data of 
N = 34 extrapolated from the region of time 200 — 2000, i.e. the linear fits of the previous 
figures. The cross at T = 1.8 is the value extrapolated from the region of time 200 — 6400 at 
N = 66. 
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4.1 The correlation functions 



Up to now we have considered equal time correlation functions. Let us see what happen to 
the correlation functions at different times. We will study the correlation C(t,t w ). 

We introduce the variable s = t/t w . According to simple aging the correlation functions 
should become a function of only s in the limit of large times. Of course the FDT region, 
which is located at finite t also when t w goes to infinity, is squeezed at s = 0, so that we 
expect that the function C becomes discontinuous at s = 0. Moreover the limit s — > give us 
information on the value of C(oo) 

limC(t,t w ) = C(oo), (34) 

where it is understood that the limit is done always remaining in the region where t » 1. 

In fig (0) we plot the correlation function at T = 1.8 for N = 66 at different values of t w 
(512, 2048, 8912) as function of s 5 = (t/t w )- 5 . 

We have plotted the data as function of s 5 and not of s for two reason: 

• To decompress FDT region at s=0 in order to see better the building of the discontinuity 
at s = 

• To expose the apparent stretched exponential decay of the correlation function for large 
values of s. 

As we can see the two regions FDT and aging are quite clear. It is also evident that C(oo) 
is different from zero at this temperature (it is obviously zero in the high temperature phase; 
we have checked this result, but for reasons of space we do not show the corresponding data). 
There is an overall drift of the correlation as function t w which seems to disappear at large 
times. This is not a surprise because we have seen that also the correlation function at s = 
show a residual dependence on t w . 

In order to see how simple aging is satisfied for other quantities we introduce the quantity 
Q(t,t w ) (introduced in ||), defined as: 

Q{t,t w )= e /(x ' (t+ ^ 2 ~ Xfc(U) (35) 

where we have chosen the function / in an appropriate way, i.e. 

'M - A- (36) 

with d = .3. The function / is very small when x » .3 and near to 1 for x < .3. 

The value of Q will be a number very near to 1 for similar configurations (in which the 
particles have moved of less than a) and it will be much smaller value (less than .1) for 
unrelated configurations; using the same terminology as in spin glasses |16|, £|, |5|] Q can be 
called the overlap of the two configurations: with good approximation Q counts the fraction 
of particles which have moved less that d. 

The data for Q are shown in fig. We see that also in this case we have some violation 
of simple scaling, but the violations are definitely smaller, also because, due the its definition 
the quantity Q is normalized to 1 at s = 0. 
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Figure 7: On a logarithmic scale we see the correlation function at F = 1.8 for N = 66 at 
different values of t w (512, 2048, 8912) as function of s 5 = (t/t w )- 5 . 
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4.2 The response function 



We have followed a standard procedure [17, 18| to measure the off-equilibrium response func- 



tion in simulations: we have kept the system in presence of an external field e up to time t w 
and we have removed the field just at this time. 

If e is sufficient small, we have that the stress as function of time, is related to the integrated 
response R by the relation 

a ( t+ J w) = S(t + t w ) = R(t + t w , 0) - R(t, t w ) (37) 

As far as the limit of R(t, t w ) when t — > oo does not depend on t w the quantity defined in eq. 
(|37D goes to zero when t — > oo. In this way we can get the value of the integrated response by 
measuring the stress density as function of time 

In our case (where we use the stress as a perturbation) the physical interpretation of the 
procedure is quite clear. We start by putting the systems in a box which is not cubic (because 
e / 0), but two sides are slightly longer of the third. At time t w we change the form of the 
box to a cubic one. In this way we deform the the system and we induce a stress which will 
be eventually decay. In the high temperature phase, where the system is liquid, the stress will 
disappear in a short time. On the contrary, in the glassy phase, we shall see that the stress 
remains for a much longer time (as expected ||) and it shows an interesting aging behaviour. 

The choice of e is crucial. In principle its value should be infinitesimal. However the signal 
is proportional to e while the errors are e independent. The errors on the response function 
grow as e -1 . On the other hand if we take a too large value of e we enter in a non linear 
region. The analysis of the non linear dependence of the stress as function of e would be very 
interesting, but it goes beyond the scope of this paper. Here we restrict to the linear region. 
We have taken data at e = .1 and e = .05 and we have seen that there are some non-linear 
effects. No non-linear effects have been detected at e = .02 and e = .01. All the data we 
present in this paper come from e = .01 and they are reasonable free of systematic effects. 
As a further check we have compared the value of Q measured at e = 0, as in the previous 
subsection, and the value of Q at e = .01 and they differs by less than 1%. 

In fig. (||) we plot the response function S at T = 1.8 for N = 66 at different values of t w 
(128, 512, 2048) as function of s' 5 = (t/t w )' 5 . In this case aging show much smaller violations 
than in the case of the correlations (due the fact that the value at s = is much less dependent 
on time). Also in this case we see that the building of a discontinuity at s = followed by a 
smooth function of s. 



4.3 Fluctuations and response together 

The crucial step would be now to plot the fluctuations and the response together. The previous 
equation tell us that 

W = X(C) (38) 

so that it is convenient to plot S versus C. The slope of the function S(C) is thus X(C). 
In the one step replica symmetry breaking scheme we expect that: 

S = mC for C<C(oo), 
S = C + (1 - m)C(oo) for C<C(oo). (39) 
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Figure 9: On a logarithmic scale we see the response function at T = 1.8 for N = 66 at 
different values of t w (128,512, 2048, 8192) as function of s- 5 = (t/t w \ 5 . 
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Figure 10: The response S as function of C at T = 1.8 for N = 34 at different values of t v 
(128, 512, 2048, 8912). 
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In fig. flTnp we find the data for the response S as function of C at T = 1.8 for iV = 34 
at different values of t w (128, 512, 2048, 8912). We can see that also here there is strong 
dependance on t w . We can distinguish two regions: 

• A high C region (typically C > C(oo)), where the function 5* has nearly slope one. 

• A region of smaller C, where the dependence of S on C is more complex. 

We remark that the data for smaller value of the waiting time seems to be very similar to the 
predictions of the one step replica theory but this conclusion is not so clear for larger values 
of the waiting time. 

In fig. (II) we find the data for the response S as function of C at T = 1.8 for N = 66 
at different values of t w (128, 512, 2048, 8192). We find a behaviour qualitatively similar 
to the one with a smaller number of particles at small waiting time; however here the data 
at t w = 2048 and t w = 8192 are much more similar, giving us the hope to be near to the 
asymptotic limit in the case of this larger system. 
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Figure 11: The response S as function of C at T = 1.8 for iV = 66 at different values of t. 
(128 512, 2048, 8192). 
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In order to expose the existence of a region where the FTD is valid |18[] , it is convenient to 
define the function D = C — S . It is evident that in the FDT region the function D(C) must 
be equal to a constant. In order to test this prediction and to verify the existence of a region 
where the FDT relation holds, in fig. (O) we plot the the function D versus S at V = 1.8 for 
iV = 66 at different values of t w (512, 2048, 8912). A plateaux region is quite evident. The 
level of the plateaux is still dependent on t w . The straight line correspond to one step replica 
symmetry breaking with m — .3. 

The behaviour of the function S(C) at small C would be quite interesting to assess. We 
have three possibilities of increasing complexity 
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Figure 12: The difference D = C — S as function of C at V — 1.8 for N = 34 at different values 
of t w (512, 2048, 8912). The straight line correspond to one step replica symmetry breaking 
with m = .3. 

• The function S(C) is zero for C < C(oo). This correspond to the simple situation where 
replica symmetry is not broken. 

• The function S(C) is linear for C < C(oo) with slope m. This correspond to one step 
replica symmetry breaking. 

• The function S(C) has a power behaviour for C < C(oo) with exponent greater than 
one. In this case the replica symmetry is broken in a continuous way. 

The previous data for S(C) plotted in double logarithmic scale may be useful to clarify 
this point, see fig. QTB"!). We find that the data for small t w seem to be linear in agreement with 
one step replica symmetry breaking predictions. There is a rather strong dependence on t w in 
this region and the extrapolation to t w = oo cannot safely done. A value of m around .3 (or 
more) is compatible with the data. This is an interesting problem that must be investigated 
further. 

5 Conclusions 

It has been shown in the previous section that there is a low temperature region where the 
large time extrapolation of the energy and of the stress autocorrelation function can be done by 
assuming simple power low corrections. The extrapolated value quite likely do not correspond 
to equilibrium values, but to metastable values and the real equilibrium values are reached 
only at much bigger time. 

The behaviour in this region is quite different from the high temperature region. The 
autocorrelation function of the stress is bigger than the analytic continuation of its value 
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from the high temperature region. Approximate simple aging is found. The function X(C) 
of the CK theory has been computed. Also this function shows a dependence on the value 
of t w which seems to decrease by increasing t w . The extrapolation of the function X(C) at 
infinite (or very large) t w it is a delicate problem, that should be treated in a careful way. 
At first sight is seems likely that in order to get a reliable extrapolation one need to increase 
the time t w by a factor 10, increase somewhat the statistics in order to decrease the errors 
and simultaneously go to larger samples in order to avoid finite volume effects, which become 
stronger and stronger when the time increase []. One probably needs a factor 100 more in 
computer time, which is not a very large amount (considered that this computation has taken 
a few weeks of a workstation), but it would certain go outside the exploratory aim of this work 
and of the capabilities of the hardware we have used. 




Figure 13: The response S as function of C at T = 1.8 for iV = 66 at different values of t w 
(128 512, 2048, 8192) on a double logarithmic scale. The straight line correspond to a linear 
behaviour. 

It is likely that something smarter can be done. Two possibilities are immediate: 

• We can repeat the same procedure where the temperature is slowed cooled in steps of 
total length t w starting from a finite temperature to the final temperature. In this way 
is possible that aging and the other phenomena survive, but the finite time corrections 
could be much smaller JL9]. 

• We could change the quantity we study. It may be possible that the stress is not the 
best suited quantity to investigate; indeed the quantity Q seems to satisfy much better 
simple aging with much smaller corrections. 

1 We should also increase the value of s in order to explore better the region where C becomes small. 
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It seems to me that there is a quite large amount of phenomena to investigate in the 
approach to equilibrium of glasses. It would be interesting to see if the indications given in 
this paper will be confirmed by more accurate and lengthy computations. 
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Appendix: spin glasses 

In this appendix I will briefly recall the situation for spin glasses. Aging and the predictions 
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Figure 14: The value of q as function of the temperature in a 4-dimensional spin glass. 



of the CK theory have been carefully analyzed in [|l!J: a comparison of their results with our 
would be instructive. Here we limit ourself to the following general remark. 
In spin glasses the relevant quantity is the total magnetization 

M = J>. (40) 

i 

It is easy to see that at zero magnetic field in the random bond model (due to gauge invariance) 
the following relation holds at all temperatures. 

< M W >' = 1. (41) 
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We can now compute in simulations the magnetic susceptibility (x), i.e. by inserting 
a small magnetic field h and my measuring the ratio < M > jh. Following the previous 
discussion (^) we have that 

X = /?(!-?)■ (42) 
An example of the function q computed in simulations for spin glasses (a four dimensional 



model) is showed in fig. fll4|) , where the data are taken from |19|| . 

There is a striking difference from spin glasses and the case of binary glasses presented 
here, which it is worthwhile to note: 

• In glasses the response to the stress tensor is computed from symmetry arguments and 
the transition is present in the fluctuations. 

• In glasses the fluctuations of the magnetization are computed from symmetry arguments 
and the transition is present response, (i.e. the susceptibility). 

However it seems to me that this difference does not seem to have deep physical consequence 
and it is just an effect of the different choice of observables. 
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